VII. BOUNDARY LAYER FLOWS 
The previous chapter considered only viscous internal flows. 


Viscous internal flows have the following major boundary layer characteristics: 


* An entrance region where the boundary layer grows and dP/dx # constant, 
* A fully developed region where: 


e The boundary layer fills the entire flow area. 


e The velocity profiles, pressure gradient, and Ty are constant; 
i.e. they are not equal to f(x), 


e The flow is either laminar or turbulent over the entire length of the flow, 
i.e. transition from laminar to turbulent is not considered. 


However, viscous flow boundary layer characteristics for external flows are 
significantly different as shown below for flow over a flat plate: 
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Fig. 7.1 Schematic of boundary layer flow over a flat plate 


For these conditions, we note the following characteristics: 


e The boundary layer thickness 6 grows continuously from the start of the 
fluid-surface contact, e.g. the leading edge. It is a function of x, not a 
constant. 


e Velocity profiles and shear stress t are f(x,y). 
e The flow will generally be laminar starting from x = 0. 


e The flow will undergo laminar-to-turbulent transition if the streamwise 
dimension is greater than a distance Xçr corresponding to the location of 
the transition Reynolds number Recr. 


e Outside of the boundary layer region, free stream conditions exist where 
velocity gradients and therefore viscous effects are typically negligible. 
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As it was for internal flows, the most important fluid flow parameter is the 
local Reynolds number defined as 








Res pU,.x _ U X 
u v 
where 
p = fluid density u = fluid dynamic viscosity 


v = fluid kinematic viscosity U — = characteristic flow velocity 


x = characteristic flow dimension 


It should be noted at this point that all external flow applications will not use a 
distance from the leading edge x as the characteristic flow dimension. For 
example, for flow over a cylinder, the diameter will be used as the characteristic 
dimension for the Reynolds number. 


Transition from laminar to turbulent flow typically occurs at the local transition 
Reynolds number, which for flat plate flows can be in the range of 


500,000< Re, <3,000,00 


With Xer = the value of X where transition from laminar to turbulent flow occurs, 
the typical value used for steady, incompressible flow over a flat plate is 


Ree = = 500,000 


cr 


Thus for flat plate flows for which: 
X < Xer the flow is laminar 


X > Xer the flow is turbulent 


The solution to boundary layer flows is obtained from the reduced “Navier — 
Stokes” equations, i.e., Navier-Stokes equations for which boundary layer 
assumptions and approximations have been applied. 
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Flat Plate Boundary Layer Theory 


Laminar Flow Analysis 





For steady, incompressible flow over a flat plate, the laminar boundary layer 
equations are: 


Conservation of mass: du + av = (0 
Ox dy 

Ou Ou 1l dp 19 Ou 
ry : V =.= SE So n 
X' momentum u T + Jy r et ot) 
KA , Op _ 
Y' momentum: —iE=() 

dy 


The solution to these equations was obtained in 1908 by Blasius, a student of 
Prandtl's. He showed that the solution to the velocity profile, shown in the table 
below, could be obtained as a function of a single, non-dimensional variable n 
defined as 


Table 7.1 the Blasius Velocity Profile 
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with the resulting ordinary a4 0.19894 a4 aiT 
i , as 0.26471 36 192333 
differential equation: La 0.32979 1s PETTE, 
1 12 02047A 4.0 99552 
14 ‘asa 4 (hi 
f” +=ff'=0 Lä 1.51674 44 97587 
2 E 0ST 46 0.98269 
Lal DATT 44. 28779 
u 22 0.68132 5.0 mans 
i s= 24 0.72899 : OOOO 

and T (n)= U p DTTA j 






co 


Boundary conditions for the differential equation are expressed as follows: 
aty=0, v=0 — f(0)=0 ; ycomponent of velocity is zero at y = 0 


aty=0,u=0 >f’(0)=0; x component of velocity is zero at y = 0 
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The key result of this solution is written as follows: 





With this result and the definition of the boundary layer thickness, the following 
key results are obtained for the laminar flat plate boundary layer: 











Local boundary layer thickness 6 (x) _ _5X 
V Rex 
Local skin friction coefficient: ER 0.664 
(defined below) RT Re, 
Total drag coefficient for length L ( integration Ge 1.328 
of Tw dA over the length of the plate, per unit D~ R 
€ 
area, divided by 0.5 p U») Vx 
T,, (x) F,/A 








where by definition Ce = and Cp = 


1 9) 1 2 
zP Us 3P Ua 


With these results, we can determine local boundary layer thickness, local wall 
shear stress, and total drag force for laminar flow over a flat plate. 

Example: 

Air flows over a sharp edged flat plate with L = 1 m, a width of 3 m and 
U» =2 m/s. For one side of the plate, find: (L), Cf (L), Ty(L), Cp, and Fp. 


Air: p = 1.23 kg/m” v=1.46 E-5 m’/s 
U „L 2m/s*215 
First check Re: Re, = —2~ =< "> = 294, 520 < 500,000 


Key Point: Therefore, the flow is laminar over the entire length of the plate and 
calculations made for any x position from 0 - 1 m must be made using laminar 
flow equations. 
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Boundary layer thickness at x = L: 

SL _ 5*2.15m 
Re, 4294, 520 
Local skin friction coefficient at x = L: 


0.664 0.664 


Re, 294,520 


Surface shear stress at x = L: 





ô(L)= = 0.0198 m = 1.98 cm 


= 0.00122 





C(L) = 


T, =1/2 pU} C, =0.5*1.23kg/m *2° m’ /s° *0.00122 
T,, = 0.0030 N/m’ (Pa) 


Drag coefficient over total plate, 0 — L: 
1.328 E 1.328 


Re, (294,520 


= 0.00245 





C (L) = 


Drag force over plate, 0 — L: 


F, =1/2pU} C, A=0.5*1.23kg / m *2° m /s° *0.00245*2*2.15m° 
F, =0.0259N 


Two key points regarding this analysis: 
1. Each of these calculations can be made for any other location on the plate 
by simply using the appropriate x location for any x < L. 


2. Be careful not to confuse the calculation for Cr and Cp. 


Cr is a local calculation at a particular x location (including x = L) and 
can only be used to calculate local shear stress at a specific x, not drag 
force. 


Cp is an integrated average over a specified length (including any x < L) 
and can only be used to calculate the average shear stress over the entire 
plate and the integrated force over the total length. 
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Turbulent Flow Equations 





While the previous analysis provides an excellent representation of laminar, flat 
plate boundary layer flow, a similar analytical solution is not available for 
turbulent flow due to the complex nature of the turbulent flow structure. 


However, experimental results are available to provide equations for key flow 
field parameters. 


A summary of the results for boundary layer thickness and local and average skin 
friction coefficient for a laminar flat plate and a comparison with experimental 
results for a smooth, turbulent flat plate are shown below. 


Laminar Turbulent 


0.16x 
oe RE ôl) = oT 











0.664 c, 0027 
A E 
C= 1.328 c 0.031 for turbulent flow over 
D~ = 7 tire plate, 0 — L, i.e. 
Re D 1/7 enure plate, , 
Re; assumes turbulent flow 
in the laminar region 
or 
C ei y? local drag coefficient based on local 
2 Q9. wall shear stress (laminar or turbulent 
ee pP ll sh (lami bul 
flow region) 
and 





Cp = total drag coefficient F/A z ( > )' L 
based on the integrated Cp ae! p U2 = (PUn A l Tw (x)wdx 
2, co 


force over the length 0 to L 
A careful study of these results will show that, in general, boundary layer 


thickness grows faster for turbulent flow, and wall shear and total friction drag are 
greater for turbulent flow than for laminar flow given the same Reynolds number. 
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It is noted that the expressions for turbulent flow are valid only for a flat plate with 
a smooth surface. Expressions including the effects of surface roughness are 
available in the text. 


Combined Laminar and Turbulent Flow 
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Flat plate with both laminar and turbulent flow sections 


For conditions (as shown above) where the length of the plate is sufficiently long 
that we have both laminar and turbulent sections: 


* Local values for boundary layer thickness and wall shear stress for either 


the laminar or turbulent sections are obtained from the expressions for (x) 
and Cf, for laminar or turbulent flow as appropriate for the given region. 


* The result for average drag coefficient C p and thus total frictional force 


over the combined laminar and turbulent portions of the plate is given by 
(assuming a transition Re of 500,000) 


0.031 0.031 
Rey’ (5x10°) 





Cp = 


* Calculations assuming only turbulent flow can typically be made for two cases: 


1. when some physical situation (a trip wire) has caused the flow to be 
turbulent from the leading edge or 

2. ifthe total length L ofthe plate is much greater than the length Xer of 
the laminar section such that the total length of plate can be considered 


turbulent from x = 0 to L. Note that this will over predict the friction 
drag force since turbulent drag is greater than laminar. 


With these results, a detailed analysis can be obtained for laminar and/or turbulent 
flow over flat plates and surfaces that can be approximated as a flat plate. 
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Figure 7.6 in the text shows results for laminar, turbulent and transition regimes. 
Equations 7.48a & b can be used to calculate skin friction and drag results for the 
fully rough regime. 


-2.5 
Cr = (287 +1.581ogž ) (7.48a) 


L -2.5 
Cp = (1.89 + 1.62 log J (7.48b) 


Equations 7.49a & b can be used to calculate total Cp for combined laminar and 
turbulent flow for transition Reynolds numbers of 5x10 and 3x10° respectively. 








1 1440 
D% 1 oe E Re € trans 7 5010" 
Re; Re, 
0.031 8700 
b= >r Rean, = 3x10° 
Re; Re, 


Example: 


Water flows over a sharp flat plate 2.55 m long, 1 m wide, with U.. = 2 m/s. 
Estimate the error in Fp if it is assumed that the entire plate is turbulent. 


Water: p = 1000 kg/m” v = 1.02 E- m/s 


* 
Reynolds number: Re, = guh aooo 5E 6 > 500, 000 
V mIs 


with Re, = 500,000 = x,, =0.255m (or 10% laminar) 


a. Assume that the entire plate is turbulent 


0.031 0.031 
Cos- ma = ea 0008423 
Re; (5x ° 
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kg 2m 


Fp = 0.5pU,, Cp A =0.5. WOOO ae 2 -0.003423-2.55m 


F- =17.46N This should be high since we have assumed that 
D i the entire plate is turbulent and the first 10% is 
actually laminar. 


b. Now consider the actual combined laminar and turbulent flow: 


_ 0.031 1440 0.031 1440 


= Sa a 0003135 
Re!” Re, (x10°) 5x10° 


D% 


Note that the Cp has decreased when both the laminar and turbulent sections are 
considered. 


k 
Fp =0.5PU2 Cp A=0.5-1000-5 S 2a -0.003135 -2.55 m? 
Fp =15.99N { Lower than the fully turbulent value} 
17.46-15. 
Error = U -100 = 9.2% high 
15.99 


Thus, the effect of neglecting the laminar region and assuming the entire plate is 
turbulent is as expected. 


Question: Since Xer = 0.255 m, what would your answers represent if you had 
calculated the Re, Cp, and Fp using x = Xer = 0.255 m? 


Answer: You would have the value of the transition Reynolds number and the 
drag coefficient and drag force over the laminar portion of the plate 
(assuming you used laminar equations). 


If you had used turbulent equations you would have red marks on your 
paper. 
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Von Karman Integral Momentum Analysis 


While the previous results provide an excellent basis for the analysis of flat plate 
flows, complex geometries and boundary conditions make analytical solutions to 
most problems difficult. 


An alternative procedure provides the basis for an approximate solution which in 
many cases can provide excellent results. 


The key to practical results is to use a reasonable approximation to the boundary 
layer profile, u(x,y). This is used to obtain the following: 


ô 
a. Boundary layer mass flow: m= lp ubdy 
0 


where b is the width of the area for which the flow rate is being obtained. 


du 
b. Wall shear stress: T, = LU 
dy = 


You will also need the streamwise pressure gradient da for many problems. 
X 


The Von Karman integral momentum theory provides the basis for such an 
approximate analysis. The following summarizes this theory. 


Displacement thickness: 





Consider the problem 
indicated in the adjacent 
figure: y 





A uniform flow field with f Streamline _ 






velocity U. approaches a 
solid surface. As a result 

of viscous shear, a ; 
boundary layer velocity Simed t 


profile develops. "| : effect 


| 
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A viscous boundary layer is created when the flow comes in contact with the solid 
surface. 


Key Point: Compared to the uniform velocity profile approaching the solid 
surface, the effect of the viscous boundary layer is to displace 
streamlines of the flow outside the boundary layer away from the 
wall. 


With this concept, we define § = displacement thickness 


x 
6 = distance the solid surface would have to be displaced to maintain the 


same mass flow rate as for non-viscous flow. 


From the development in the text, we obtain 


Thus, the displacement thickness varies only with the local non-dimensional 


* 
velocity profile. Therefore, with an expression for u/ Us , we can obtain ò = 


(8). 


Example: 


2 
Given: a = z) — (2) determine an expression for 5 = f(d) 
U ô ô 


Note that for this assumed form for the velocity profile: 


co 


1.At y =0, u=0 © correct for no slip condition 


2.At y = 6, u = Us correct for edge of boundary layer 
3. The form is quadratic 


To simplify the mathematics, 
let n=y/d, aty=0, n=0; aty=d, n=1; dy=ddn 


u 2 
Therefore: —— =2n-— 
erefore U n-n 


co 
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~ 1 2 3)1 
Substituting: ô =|Q-2n+n )Sdn=6 p 
i 0 


* o] 
which yields ô => 


Therefore, for flows for which the assumed quadratic equation approximates the 
velocity profile, streamlines outside of the boundary layer are displaced 
approximately according to the equation 


This closely approximates flow for a flat plate. 


Key Point: When assuming a form for a velocity profile to use in the Von 
Karman analysis, make sure that the resulting equation satisfies both surface and 
free stream boundary conditions as well as has a form that approximates u(y). 


Momentum Thickness: 





The second concept used in the Von Karman momentum analysis is that of 
momentum thickness - 0 
The concept is similar to that of displacement thickness in that O is related to the 


loss of momentum due to viscous effects in the boundary layer. 


Consider the viscous flow 


regions shown in the uf A ad Di A A= z 


adjacent figure. 






2 


Define a control volume as r 
D (on fluid) 


shown and integrate around 





B 


: 
the control volume to obtain / Simulated i 
the net change in |. , ase 
momentum for the control 

volume. 
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If D = drag force on the plate due to viscous flow, taking the fluid as the 
control volume, we can write 


-D = ) (momentum leaving c.v. )- }, (momentum entering c.v. ) 


Completing an analysis shown in the text, we obtain 


ô u u 
P l-ie 
D/A 
Using a drag coefficient defined as Co =] 2 
2P Ua 
20(L 
We can also show that Cp = < ) 


where: @(L) is the momentum thickness evaluated over the length L. 


Thus, knowledge of the boundary layer velocity distribution u = f(y) also allows 
the drag coefficient to be determined. 


Momentum integral: 





The final step in the Von Karman theory applies the previous control volume 
analysis to a differential length of surface. Performing an analysis similar to the 
previous analysis for drag D we obtain 








This is the momentum integral 


Ty m 5 U dU d VU: o) for 2-D, incompressible flow 
dx 





p “y + and is valid for laminar or 
turbulent flow. 
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z dU ô dP Therefore, this analysis also accounts for 
where 0 U,,—= =————___ the effect of freestream pressure 
dx p dx gradient. 


For a flat plate with non-accelerating 
flow, we can show that 





= const., 7 BE 
x 





P = const., U 


co 


Therefore, for a flat plate, non-accelerating flow, the Von Karman momentum 
integral becomes 


From the previous analysis and the assumed velocity distribution of 


2 
BONN ape 
g745)-(5) = 20-7 


The wall shear stress can be expressed as 


d 2 on 2 uU 
7 dy w y=0 


Also, with the assumed velocity profile, the momentum thickness © can be 
evaluated as 


ô 
ree pee pees dy 
oU U 


or 
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6 ô 
o= @n-n*)l -21 +n Ban= = 


We can now write from the previous equation for Ty 


wat So ng Pe 
Ee fe 15" dr 


Equating this result to Eqn. A we obtain 
2 dô 2uU 
pa pi e 
15 dx Ô 





or 
15 
ddd= H dx which after integration yields 
pu,, 
5 _ (30m -= ya S4 
= ———_. r = 
pU, ° Re, 


Note that this result is within 10% of the exact result from Blasius flat plate 
theory. 


Since for a flat plate, we only need to consider friction drag (not pressure drag), 
we can write 
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Exact theory has a numerical constant of 0.664 compared with 0.73 for the Von 
Karman integral analysis. 


It is seen that the Von Karman integral theory provides the means to determine 
approximate expressions for 


6, Ty, and Cr 


using only an assumed velocity profile. 


Solution summary: 


. Assume an analytical expression for the velocity profile for 
the problem. 


. Use the assumed velocity profile to determine the solution 
for the displacement thickness for the problem. 


. Use the assumed velocity profile to determine the solution 
for the momentum thickness for the problem. 


. Use the previous results and the Von Karman integral 
momentum equation to determine the solution for the 
drag/wall shear for the problem. 





and 7.3 and does not have to be accounted for separately. 
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